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ABSTRACT 

Recently, detections of a high-energy 7-ray source at the position of the Galactic center have been 
reported by multiple gamma-ray telescopes, spanning the energy range between 100 MeV and 100 
TeV. Analysis of these signals strongly suggests the TeV emission to have a morphology consistent 
with a point source up to the angular resolution of the H.E.S.S. telescope (approximately 3 pc), while 
the point-source nature of the GeV emission is currently unsettled, with indications that it may be 
spatially extended. In the case that the emission is hadronic and in a steady state, we show that 
the expected 7-ray morphology is dominated by the distribution of target gas, rather than by details 
of cosmic-ray injection and propagation. Specifically, we expect a significant portion of hadronic 
emission to coincide with the position of the circum-nuclear ring, which resides between 1-3 pc from 
the Galactic center. We note that the upcoming Cherenkov Telescope Array (CTA) will be able to 
observe conclusive correlations between the morphology of the TeV 7-ray source and the observed gas 
density, convincingly confirming or ruling out a hadronic origin for the 7-ray emission. 

Subject headings: (ISM:) cosmic rays — gamma rays: theory — gamma rays: observations 



1. INTRODUCTION 

Since COS-B and EGRET first observed a bright 7- 
ray source spatially coincid e nt with the Galactic cen- 
ter ( GC) ([Blitz et all 119851: iMaver-Hasselwander et~aTl 
I1998T ). subsequent observations by the Large Area Tele- 
scope onboard the Fermi space observatory (Fermi-LAT) 
and by ground-based (imaging Atmospheric Cherenkov 
Telescope, ACT) 7-ray telescopes have repeatedly ob- 
served the GC region at energies spanning 100 MeV to 
100 TeV. Unlike both radio and X-Ray observations of 
the GC, no variability has been observed in the high- 
energy regime, potentially indicating that the 7-ray emis- 
sion mechanism differs substantially f rom the low en- 
ergy regime (jChernvakova et all 120111 : lAharonian et al.l 
2009). This distinction is especially stark in light of 
the order of magnitude increase in X-ray activity from 
the region observed in 2005 , whic h was undetected in 
7-ray data ([Aharonian et al.ll2008l ). This might be ex- 
plained by models which generate the low-frequency (X- 
ray, IR, radio) emission very close to the central black 
hole (BH), while 7 rays are produced farther away from 
the BH by high-energy protons inelastically scattering 
off of the interstellar medium, a framework originally 
posited by lAharonian fc Neronovl {2005), and later by 
iLiuet alJ (|2006allU >. 

The strongest limits on the morphology of the gamma- 
ray so urce HESS J1745-290 are provided bv lAcero et al.1 
(2010) which used optical cameras mounted on each 
H.E.S.S. dish to calibrate the high-energy array, reducing 
the systematic error in the s patial pointing of th e tele- 
scope. Using this technique, lAcero et al.l (J20101) found 
the source HESS J1745-290 to be spatially coincident 
with the radio source Sgr A*, within a total error ra- 
dius of only 13y excluding the supernova remnant Sgr A 
East as the dominant source of 7-ray emission. However, 



the pulsar wind nebula G359.95-0.04 cannot be ruled out 
as a possible source of the high energy emission, and it 
may play a secondary or even a dominant role in the 
7-ray emission observed from the galactic center. We 
note, however, that this possibility would not affect the 
analysis and conclusions presented here if the gamma-ray 
emission from the PSR was produced via the injection 
of high-energy protons into the galactic medium. HESS 
J1745-290 is found to be consistent with a point source to 
within approximately 1.2' at the 95% confidence leveQ. 
While this point source resides within a diffuse 7-ray 
background which is itself centered around Sgr A*, the 
total residual flux from the inner 0.1° is constrained to be 
less th an 15% of t he point-source fl u x (lAharonian et al.1 
I2006a|) . Work by lAharonian et all (|2009f ) further con- 
strained this diffuse emission to be relatively indepen- 
dent of the GC distance - indicating it may stem from 
cosmic-ray background events rather than a diffuse signal 
corresponding to the GC. 

The LAT has also observed a point source, 1FGL 
J174 5.6-2900c, spat ially coincident with the GC re- 
gion (|Abdo et a l. 2010), although a conclusive identifica- 
tion with Sgr A* was impossible due to the relatively low 
angular resolution of the instrument (~ 0.2° at 10 GeV). 
There are indications, however, that the GC signal ob- 
served by the Fermi-LAT is difficult to model as a simple 
point so urce. While ind e pende nt po int-source analyses 
by bo th iBovarskv et al.l (|2011f ) and Chcrnvakov a et al.l 
(|2011h produced best-fit models assuming a point - source 
emission, the analysis of iHooper fc Goodenoughl (|2011f ) 
found that the 7-ray emission extends spherically out 
to approximately 50 pc from the position of Sgr A*, 
falling off with a power- law of approximately r~ 26 . 

1 This corresponds to 2.96 pc under the assumption of a solar 
position 8.5 kpc from the GC. 



iHooper fc LindeiJ (|2011l) pointed out two potential issues 
with GC point-source models at Fermi energies. First, 
a point-source best fit systematically over-estimates the 
source flux in the presence of an appreciable diffuse back- 
ground. Secondly, even after accounting for surround- 
ing point-sources and Galactic plane emission, the ex- 
cess emission observed in the GC and not associated to 
the (best-fitted) point source, exceeds the point source 
luminosity by approximately a factor of three - in stark 
contrast to H.E.S.S observations at TeV energies. This 
second finding is con sistent with the observation by 
iBovarskv et alJ (|2011f l that the log-likelihood of the fit 
increased by 25 with the addition of a spherical symmet- 
ric parameter describing the spatial extension. 

In addition to these morphological inconsistencies, the 
origin of the GC source spectrum at high-energy 7-ray 
frequencies is puzzling. A relatively hard spectrum in the 
O.I-I GeV range significantly softens in the 1-100 GeV re- 
gion, then hardens again at TeV energies, before cutting 
off above 10 TeV. Leptonic models h ave been proposed 
to explain both the TeV emiss ion (|Atovan fe Dermen 
l2004t Irlinton fc Aharonianl 120 07) and the GeV emission 
(Kusunosc & Takahara 2012). However, no theoretically 
compelling scenarios have been proposed to explain the 
entire gamma-ray spectrum a single leptonic source pop- 
ulation^ 

Ballantv ne et al.1 (|2007T> examined the possibility that 
the TeV signal observed by HESS could be explained 
by a significant emission of hard protons from Sgr A*, 
producing the gamma-ray emission through interactions 
with gas in the galactic center region, s imilar to the ear- 
lier models of Aharoni an fc Neronovl (|2005|) . In order 
to examine the possibility that this emission be con- 
fined to the region near Sgr A*, they employed mod- 
els for the morphology of Hydrogen gas in the GC from 
iRockefeller et a l. (2004), which include a large "ring" of 
overdense gas surrounding the GC at distances from 1- 
3 pc. This dense structure, aptly named the circumnu- 
clear ring (hereafter CNR), has long been observed by 
far-infrared instruments as a torus-shaped structure in- 
clined 20° with respect to the Galactic plane and sur- 
rounding a relative underdensity of gas within the in - 
ner pc around the GC (see e.g. iBecklin et al.l I1982T ). 
iBallantvne et al.l (|2007|) then calculated the diffusion 
of charged protons within tangled magnetic fields with 
strength proportional to the local Hydrogen density. Us- 
ing these models, they calculated the expected interac- 
tion probability between high energy protons and the 
CNR as a function of energy, finding approximately 73% 
of protons emitted from the GC to encounter the CNR for 
protons with energies between 1-2.5 TeV. However they 
found this interaction probability to fall quickly with in- 
creasing proton energy due to the increasing gyroradius 
of high energy protons - only 47% and 5% of protons 
are found to encounter the CNR at energies of 10 and 
100 TeV respectively. Due to the extreme oyerden sity 
of hydrogen gas in the CNR, IBallantvne et al.l (|2007l ) as- 
sumed that protons encountering the CNR would lose en- 
ergy quickly through pion collisions, while protons which 
avoid the ring would lose energy at larger radii, produc- 
ing an exte nded gamma-ray emission at higher energy. 

Recently iChernvakova et all (|2011f) further examined 
the hadronic scenario, additionally positing that the 
whole GeV- TeV 7-ray emission may be explained by the 



injection of high energy protons at the GC. This inter- 
pretation has several natural advantages in explaining 
the entire 7-ray spectrum, including: (1) the hard spec- 
trum below 1 GeV is naturally explained by the inability 
of protons with kinetic energy below ^300 MeV to pro- 
duce 7T° in p-p collisions, (2) the bump in low energy 7- 
ray emission is produced by diffusively trapped protons 
which lose significant energy while propagating through 
the GC region, (3) the fiat spectrum at TeV energies 
is explained by protons which propagate rectilinearly 
through the GC region, without losing significant energy 
in p-p collisions, which provides a convincing match to 
the ~E -2 7-ray spectrum by using the E~ 2 proton in- 
jection spectrum derived from first-order Fermi acceler- 
ation, and (4) the turnover between the GeV and TeV 
emission is naturally explained as the transition between 
diffusive and rectilinear motion - creating a sof t spectrum 
betwe en these two energy regimes. IChernvakova et al.1 
(|2011|) also consider the possibility that two distinct pro- 
ton populations are responsible for the low- and the high- 
energy 7-ray emission. 

In the present study we focus on the morphology of 
the GC 7-ray source. Specifically, we introduce a de- 
tailed model of the interstellar gas distribution near the 
GC, and demonstrate that the 7-ray morphology from 
hadronic emission is determined primarily by the distri- 
bution of target gas rather than by the parameters de- 
scribing the diffusion of high-energy protons. We find 
that in any scenario for cosmic ray diffusion, the bulk 
of the high-energy emission falls within the point-spread 
function of all current 7-ray telescopes. However, we 
note that this will not be the case for the proposed 
Cherenkov Telescope Array (CTA), which will have the 
angular resolution required to observe a morphology in 
the 7-ray source which sha dows the observed gas density 
(|CTA Collabor ation 201l|). The outline of this study is 
as follows: Sec. [5] describes the model we employ for the 
GC gas density; Sec. Ogives qualitative arguments to de- 
scribe cosmic ray proton propagation in the GC region, 
while in Sec. [4] and [5] we present our results for the pro- 
ton population and for the 7-ray emission in hadronic 
models for the GeV- TeV emission. 

2. GAS DENSITY NEAR THE GALACTIC CENTER 

Detailed knowledge of the gas density is critical to ac- 
curately describing hadronic 7-ray emiss ion from t he GC 
region. Here we employ the model of iFerrierd (2012), 
valid for the inner ~10 pc around the GC. This model 
not only includes a diffuse halo, but also contributions 
from the CNR, Sgr A East, M-0.13-0.08 and M-0.02-0.07 
as well as bridges and streamers connecting them. We 
make several necessary simplifications to the model be- 
fore utilizing it within our numerical code. First, we 
create a spherical model by calculating the average gas 
density at each radius r. In doing this, we take a central 
value for the position of M-0.13-0.08 and make geometric 
approximations for the position of the va rious str e amers , 
keeping the total volume consistent with IFerrierd (J2012). 
We note that this will have the effect of artificially mak- 
ing the 7-ray emission spherically symmetric. However, 
in the limit that both the proton diffusion is spherically 
symmetric, and the probability of multiple interactions 
between a single cosmic ray proton and the interstellar 
medium is low, the total emission from within a given 
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Fig. 1. — Gas density (nj/cm 3 ) averaged over solid angle as a 
function of radi us from the GC (pc) obtained from a com bination 
of the work of .Ferriere et al.l 1)20071 ') and Fcrricr^ H2012T ) (see the 
text). The major feature in the gas density stems from the circum- 
nuclear ring (CNR, 1.2 pc to 3.0 pc) which contributes a gas density 
nearly two orders of magnitude larger than any other structure. 

radial bin is conserved when the gas density is smeared. 
We will discuss these assumptions in detail later, and we 
will show that they are valid for the majority of the pa- 
rameter space discussed. Secondly, we ignore the thermal 
distribution of gas, which is valid in the limit where we 
only consider collisions with highly-relativistic protons. 
To extend our simulations beyond 10 pc from the GC 
in order to capture the region relevant for the Fermi- 
L AT PSF, we a dopt the model of Eqs. 18 and 19 
of IFerriere et al.l (|2007f) and impose again spherical sym- 
metry. With these assumpti ons, we obtain gas densities 
which are lower than those in lFerrierel (|2012f ) by an order 
of magn i tude at a radius of 10 pc. Since the models of 
IFerriere! ([2012D are only suffering incompleteness at radii 
greater than 10 pc (due to unmodeled sources), we as- 
sume the gas density a t a given radius to b e t he larger 
of the values quoted by IFerriere et al.l (|2007l ) or IFerriere! 
(2012). This causes us to switch between models at a 
radius of 19.2 pc. In FigureUJwe show the resulting den- 
sity of hydrogen gas in our simulation as a function of 
the radius from the GC. We note that the morphology is 
dominated by the CNR, which provides a large boost to 
the gas density between 1-3 pc. 

3. A QUALITATIVE MODEL 

In order to understand the morphology of the hadronic 
emission, we consider propagation in four limiting 
regimes, controlled by two parameters. Protons may 
propagate either rectilinearly or diffusively through the 
GC region, and they may either undergo many, or much 
less than one, collision with the surrounding gas. While 
intermediate cases are possible, we find these four limit- 
ing cases to confine the expected proton morphology. 

We first investigate the propagation of protons in the 
rectilinear regime. This can be thought of as diffu- 
sion with a mean-free-path exceeding the confinement 
region of the simulation, yielding a diffusion constant 
D = 's, m /(6T), where l S i m is the region we are con- 



sidering, and r is the propagation time out of the re- 
gion. For our 50 pc simulation and assuming relativistic 
propagation velocities, this corresponds to a minimum D 
> 7.7 x 10 29 cm 2 s _1 for particles to propagate rectilin- 
early out of the simulation zone. We can calculate the 
probability of a proton-proton interaction between the 
cosmic ray and a target gas molecule as 
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where a pp {E) is the cross-section between a proton 
with Energy E and a cold target proton, ph is the gas 
density as a function of position, and r(t) is the radial 
position as a function of time from the injection of the 
cosmic ray until it leaves the diffusion zone. In the case 
of spherical symmetry and rectilinear propagation out 
of the center of the diffusion zone, we set r(t') = ct', 
r(0) = 0, and r(t) = r max , and then using the gas den- 
sity shown in Figure [TJ infer an interaction probability 
of ^2.5% between a given proton and the target gas be- 
fore the proton escapes the 50 pc region. In this case, 
relativistic protons lose little energy as they propagate 
through the region, and the proton spectrum mirrors the 
injected spectrum. Since the number of protons in a 
spherical shell 5r is constant for rectilinear propagation, 
the resulting 7-ray morphology mirrors the gas distribu- 
tion shown in Figure [TJ In this scenario 87% of the total 
7-ray emission is concentrated within the inner 3 pc sur- 
rounding the GC. If, as expected, this scenario describes 
the propagation of the ~10 TeV protons responsible for 
the TeV 7-ray signal, then this ratio stands in excellent 
agreement with HESS observations indicating at least 
85% of 7-ray TeV emission to be contained in the in- 
ner 3 pc. We note that so long as we assume rectilinear 
propagation, this conclusion still holds in the case of a 
true 3D model of the gas density, and the final emission 
morphology is calculable by convolving the p(r) ~ r -2 
density of pro tons with the observed gas density (e.g., 
Ferriere (2012)) and integrating over the line of sight. 

We note that for the gas densities and distance scales 
modeled here, the second possible regime, where the dif- 
fusion is rectilinear but protons undergo more than one 
collision with the target gas is insignificant, since the 
probability of having even one collision is only ~2.5%. 
Thus we exclude the second regime. 

A third regime exists when protons propagate diffu- 
sively through the region, but interact with the Galactic 
gas fewer than once. In this case, we can assume that 
the energy spectrum of protons remains constant with 
radius, and calculate the density of protons to fall off 



as p{r) 



causing the number of protons within a 



radial bin to increase linearly with radius. In this case, 
we can directly calculate the morphology of the result- 
ing 7-ray emission by using the gas density shown in 
Figure [TJ but weighting the gas density by the number 
of protons in each radial bin, which varies linearly with 
r. Weighted in this manner, we calculate the average 
ny density as 410 cm -3 , and find the minimum diffu- 
sion constant for which less than 0.1 interaction occurs 
to be 2.4 x 10 27 cm 3 s _1 . In this limit, the expected 7- 
ray morphology results from convolving the gas density 
in Figure [TJ with a proton spectrum which is constant 



in radius and has an overall density p(r) 



the 7- 



TABLE 1 



Parameter 



Chcrnyakoya^t_aL (2011) Value Value Adopted Here 



Radial Size of Simulation Region 

Duration of Proton Emission 

Proton Injection Spectrum 

Hydrogen Density during proton diffusion 

Hydrogen Density for 7-ray calculation 



3 pc 

300 - 10 4 yr 
-1.9- -2.0 
1000 cm -3 
1000 cm -3 



50 pc 

10 10 yr 

-1.9 

1000 cm -3 

See Figure 1 



List of all input parameters which differ between the assumptions of Chcrnyakova ct al. (2011) and the 
present work. 



ray morphology is therefore still dominated by the CNR, 
mirroring what was found in first regime. 

Lastly, we evaluate the regime where diffusive cosmic 
rays are expected to interact multiple times with inter- 
stellar gas before diffusing out of the GC region. In this 
case, the final cosmic ray energy spectrum is not radius 
independent, and cannot be easily computed. Further- 
more, we note that in this regime, the non-spherically 
symmetric nature of the gas becomes important in de- 
termining the steady-state cosmic-ray density, as inter- 
actions which occur while protons move through partic- 
ularly high-density regions may have an effect on the 
surrounding proton distribution. However, we note that, 
so long as the injection proton energy spectrum is suf- 
ficiently steep that partially cooled high energy protons 
are subdominant to the injected low energy population, 
this cosmic ray density is unambiguously constrained to 
fall off faster with increasing radius than in the one- 
interaction regime above. In order to provide more quan- 
titative calculations, we defer to a numerical model for 
the propagation of high energy particles throughout the 
GC region. 

We note that the r _1 dependence of cosmic-rays prop- 
agating diffusively out of the galactic center, closely mir- 
rors the brigh tness profile o f gamm a-ray emission calcu- 
lated by Chcrnyak ova et aLl (|2011D (Figure 7) for cosmic 
rays of 1-10 GeV. This is expected, as these models use 
a constant Hydrogen density, creating a one-to-one cor- 
respondence between the proton density and the bright- 
ness profile. Noti ce that the diffusion con stant at 1-10 
GeV employed by [Chcrnyak ova et al.l ()2011l ) fall between 
9.5 x 10 26 and 7.6 x 10 27 cm 2 s _1 , which lands squarely 
in the regime of diffusion with fewer than one interac- 
tion per cosmic ray (using a region with radius 3 pc and 
a gas density of 1000 cm -3 ). In the lower energy bin 
of 0.1-1 GeV, the significantly smaller diffusion constant 
leads to multiple interactions between high energy pro- 
tons and the target gas, leading to a more constrained 
distribution, as predicted in our qualitative model. 

4. PROPAGATION OF HIGH-ENERGY PROTONS 

In order to simulate the propagation of protons from 
a cent ral source, we adopt the formalism of lAloisio et al.l 
(2009), which seeks to model spherically symmetric dif- 
fusion while avoiding the issue of superluminal propa- 
gation. We note that this cosmic-ray diffusio n model is 
nearly identical to that employed bv lChernvakova et al.l 
(|2011[) and we refer the reader to that work for rele- 
vant details of our simulation parameters. In Table Q] 
we list all differences be tween our model and that of 
IChernvakova et aLl (|2011l) . Most importantly, we employ 
a much larger simulation region with a 50 pc radius, in 



order to capture the entirety of emission within the Fermi 
PSF. In order to obtain a steady state diffusion solution 
over this period, we must assume that the the emission 
continues for perio d s long er than the 10 yr assumed in 
IChernvakova et aLl (|2011 ). Since we are only concerned 
with the upper limit for the source extension in this work, 
we allow the emission to continue for 10 10 yr, but note 
that the equation becomes steady state after approxi- 
mately 10 6 yr. 

We adopt a differential proton injection spectrum 
which follows a power law ~E -1,9 with an exponential 
cutoff at 100 TeV. We normalize this injection spectrum 
in order to obtain a 7-ray int ensity matching the Ferm i 
and HESS fluxes reported by IChernvakova et al.l (|201fl) 
at a radius of 3 pc, and then extend our simulation from 
this point out to 50 pc. We find that this translates to 
an injected proton luminosity of 8.8 x 10 36 erg s _1 . 

For simplicity, we utilize an average target density 
iih = 1000 cm -3 in our calculation of the final proton 
density (but not for the final 7-ray emission, for which 
we employ the target gas density given by Figure [1}. 
This is an immaterial assumption, so long as we reside 
in a limit where the average proton undergoes much less 
than one collision with gas - an assumption which holds 
throughout the vast majority of our parameter space. 
We calculate the 7-ray emission from the steady-state 
proton density distr ibution employing the formalisms of 
IKamae et al.1 (J2006I ). We note that we consider here only 
contributions from 7-rays produced directly in the p- 
p collision, and ignore possible 7-ray contributions due 
to the inverse-Compton scattering (ICS) of photons by 
leptons produced in the same collisions. This approx- 
imation is reasonable both because ICS contributes at 
significantly lower energies than the direct 7-ray chan- 
nel (and will thus be subdominant to the larger proton 
flux at lower energies), and because the magnetic fields 
are expected to peak strongly in this region, forcing the 
ICS to be highly subdo minant to synchrotron produc- 
tion (jCrocker et al.ll20Toh . We note that variations in the 
diffusion parameters are constrained by the transition be- 
tween the GeV- TeV spectrum, and yielded no qualitative 
changes in the findings outlined below. 

5. RESULTS 

In Figure [2l we show the spectral energy distribution 
of the 7-ray signal within cones corresponding to varying 
radii. We note that more than 50% of the residual emis- 
sion is found within 3 pc at all energies, and the majority 
of this emission is created between 1-3 pc from the GC, 
when the emitted protons interact with the dense CNR. 
In Figure [3l we instead show the differential (top panel) 
and integrated (bottom) emission as a function of radius 
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Fig. 2. — Spectral energy distribution of 7-ray emission within 
radial cones of various size ranging from r < 1 pc to r < 50 pc. 
We note that the emission falls off very quickly after 3 pc, with 
almost no additional intensity for additional radial bins beyond this 
point. The overall intensity is normalized such that the emission 
within 3 pc matches the Fermi and HESS fluxes as reported by 
Chcrnvakova ct al. (2011). Due to the higher n^j in our simulation 
and longer emission period, this corresponds to a lower proton 
injected luminosity of 7.0 x 10 35 erg s _1 into protons. 



at various fixed energies. In the top panel, we show the 
total emission in radial shells of 1 pc width, while in the 
bottom panel we show the integrated emission within a 
given radius. 

We note that these features correspond closely to the 
expectations from Section [3] For instance, at energies of 
10 TeV the emission is dominated by a proton population 
of approximately 100 TeV, which propagates rectilinearly 
in our simulation. We find the 7-ray emission contained 
within 3 pc and 10 pc to be 87% and 95% respectively, 
matching the expectations from rectilinear propagation. 
At energies of 10 GeV, the 7-ray signal is dominated 
by 100 GeV electrons which propagate diffusively but 
undergo Poissonian interactions. In this case we calculate 
the 7-ray emission contained within 3 pc and 10 pc to 
be 61% and 78% respectively, closely matching an r _1 
proton density. 

Finally, in Figure [3] (bottom) we show with vertical 
lines the approximate angular resolutions of front con- 
verting events from the Fermi-LAT at 10 GeV (28 pc, 
orange solid), and at 100 Ge V (18 pc, gray dott e d), HES S 
(11 pc, gray long-dashed lAharonian et al.l (|2006bl) ). 
the HESS 95% confidence limit on th e maximum 
source extension(3 pc, gray short-dashed, lAcero et al.1 
( 
( 



2010Tl) and the anticipate d proposed resolution of CTA 



CTA Collaboration! l20TTl 2.5 pc, light blue dash-dot) 



We note that the angular resolution for the Fermi-LAT 
at 1 GeV is approximately 90 pc and is outside the plot- 
ted range. For all existing telescopes, the modeled radial 
dependence could easily be confused with a point source 
located at the GC. However, similarly to HESS, CTA will 
be able to place limits on the source extension which are 
a factor of few smaller than the quoted angular resolu- 
tion. In this case, CTA will clearly detect structure in the 
7-ray signal, including, possibly, bright emission coinci- 
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Fig. 3. — Differential (top) and integrated (bottom) radial emis- 
sion for 7-ray energies of 1 GeV (solid black), 10 GeV (dotted 
red), 100 GeV (green dashed), 1 TeV (blue long dash) and 10 TeV 
(purple dot-dashed). In the bottom panel we show vertical lines 
corresponding to the angular resolution of the Fermi-LAT at 100 
GeV (orange solid), Fermi-LAT at 10 GeV (brown dotted), HESS 
(gray long dashed), CTA (projected, blue dash-dot), as well as the 
HESS 95% confidence limit on the source extension (HESS PS, 
gray short dashed). 

dent with the CNR. We expect that this signature would 
provide extremely strong evidence in favor or against the 
hadronic nature of this emission. 

6. DISCUSSION AND CONCLUSIONS 

We have shown that hadronic emission stemming from 
the inner 3 pc of the Galaxy and following a power-law 
injection spectrum compatible with Fermi acceleration 
will naturally produce an emission spectrum comparable 
to that observed by both HESS and the Fermi-LAT. We 
showed analytically and numerically that the morphol- 
ogy of 7-ray emission is determined primarily by the gas 
morphology, and the majority of the emission falls be- 
low the PSF of all current 7-ray telescopes. While this 
is in extremely good agreement with HESS observations 
reporting 85% of the TeV emission to be c onfined within 
3 pc of the GC (|Aharonian et al.lr20 06a). the hadronic 
model is currently in some tension with Fermi-LAT ob- 
servations which imply that the GC source may be ex- 
tended. 

Specifically, our models find 71% and 86% of the 1 GeV 
7-ray emission to fall within 3 pc and 10 pc from the 
GC respectively. This result stands, at face value, in 
some tension with Ferm i-LAT observations of the GC 
(iHooper fc Lin den 2011). After background subtraction, 
iHooper fc Linden! (2011) find a residual emission which 
is not well-modeled by a point source convolved with 
the Fermi-LAT PSF. If modeled as an extended emission 
from dark matter (but independent of any dark matter 
properties), they derive a best-fit power-law fall-off to 
the emission ~ r~ 2 - 6 , which indicates that roughly 32% 
and 53% of the 7-ray signal should originate within 3 pc 
and 10 pc of the GC. 

It is important to remark that it will be difficult to 
definitively conclude that hadronic models are ruled out 
by Fermi observations, as significant background sub- 
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traction is necessary in order to determine the resid- 
ual background signal which should be observed by the 
Fermi-LAT. This includes a dominant foreground from 
emission in the galactic plane which lies along the line 
of sig ht between the s olar position and the galactic 
plane (jHooper fc L inden 2011). This intensity of residual 
emission can be further complicated by uncertainties in 
the mo deling of point source contributions at the galactic 
center (|Bovarskv et al.ll2011|) . Lastly, the 7-ray morphol- 
ogy for hadronic models is neither point-like nor power- 
law with radius, making it difficult to directly compare 
these results with previous works. Still, our results in- 
dicate that a substantial source extension in Fermi-LAT 
data is not expected from hadronic emission models. 

It may appear that the conclusion above can be cir- 
cumvented if the hadronic emission from the GC is highly 
polar, allowing protons to propagate out of the central 
3 pc while avoiding the CNR. However, this is difficult 
to reconcile with TeV observations: with purely polar in- 
jection, we find that HESS would clearly observe an ex- 
tended source above and below the Galactic plane. Our 
conclusions would also be affected if the GeV emission 
mechanism were distinct from the TeV emission - either 
due to different primary proton populations or to multi- 
ple low-energy hadronic sources contributing to the GeV 
but not to the TeV emission. 

Notice that a time dependence in the hadronic emis- 
sion could be engineered to provide an energy-dependent 
source extension. Further multiwavelength observa- 
tions will be necessary to constrain this scenario. An- 
other alternative pertains to the possibility that the 
CNR is not a stable feature of the galactic-center re- 
gion, but is instead transient in nature. While stud- 
ies of galactic tracers such as HCN and HCO + found 
extremely high gas masses of approximately 10 6 M , 
which would be stable agai nst tidal disruption by 
the supermassive-black hole ((Christopher et all 120051 : 
iMontero-Castano et al.1 I2009I) . very recent observations 
using the GR EAT telescope to p erform a CO excitation 
analysis (Re quena- Torres et al.l 2012) obtain a best fit 
which reduces the total mass of the CNR by approxi- 
mately two orders of magnitude, and thus find the CNR 
to be susceptible to tidal disruption, implying that the 
feature could be transient. We note that the density 
us ed in thi s pape r, based on the central value reported 
bv lFerriera (|2012l ) stands at 2 x 10 5 M Q , in between these 
extreme values. The implication of a transient CNR in 
the context of this work is difficult to positively deter- 
mine, as it would depend sensitively on the effect of the 
CNR not only on the gas density, but on the diffusion 
parameters of the region as well. If the latter effect is 
neglected, then it is possible that the fall-in of the CNR 
would lead to enhanced low-energy emission, as protons 
at this energy have built up without losing significant en- 
ergy until the CNR moves into regions with higher proton 
density nearer the GC. 

We note that while our work considers a similar 
model for the gas de nsity near the GC to that of 
iBallantyne et al.l (p007). we obtain an energy depen- 
dence for the gamma-ray morphology which is qua lita- 
tively different. Specifically, Ballantv ne et al.l (|2007l ) cal- 
culates the probability of a cosmic ray proton entering 
the CNR before it diffuses out of the target region. Since 
this probability will decrease as the protons mean free 



path becomes larger, they determine that the gamma- 
ray morphology should become more diffuse at high en- 
ergy. However, this implicitly assumes that gamma-rays 
moving through the CNR efficiently lose their energy be- 
fore leaving the region. However, empl oying the mode l 
for Hydrogen in the CNR determined bv lFerrierei (J2012), 
and using a model where the propagation of TeV pho- 
tons is rectilinear on the distance scale inhabited by the 
CNR, we find that protons escape from the CNR with 
only negligible energy losses. Additionally, we find that 
these energy losses remain negligible even through the 
majority of the diffusive regime. 

This creates a significant difference in the neces- 
sary proton injection spectra necessary to match the 
HESS dataset in thi s work, compared to that of 
IBallantyne et al.l ([2007) . In our work the existence of rec- 
tilinear proton propagation throughout the HESS regime 
implies that the injected E -19 proton produces an al- 
most equivalent E~ 19 gamma-ra y spec trum. However, 
in the work of IBallantyne et al.l (|2007f ) the decreasing 
interaction probability between high energy protons and 
the molecular gas requires a much harder injected proton 
spectrum (E -075 ) to account for the observed gamma- 
ray spectrum. Finally, we note that while the percent- 
age of high energy emission contained within the inner 
3 pc in our model is similar to the probability of a 1 
TeV proton encounteri ng the CNR as determined by 
IBallantyne et al.l ([2007), this is fairly coincidental: the 
relative emission in our model depends sensitively on the 
Hydrogen gas density outside of the i nner 3 pc, which 
is not taken to account in the model of IBallantyne et all 
(|200l . 

Finally, we note that the greatly improved angular 
resolution of CTA will provide a much crisper picture 
than current instruments - potentially leading to the 
observation of both a central deficit of 7 rays in the 
central pc from the GC, and of a bright ring of emis- 
sion coincident with the position of the CNR if the 
GC source is of hadronic nature. This morphology 
would be difficult to explain with leptonic or dark 
matter scenarios, which would not be expected to 
correlate with the local gas density. Thus we believe 
that future CTA observations will convincingly confirm 
or rule out the hadronic emission scenario presented here. 
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